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1. HAMILTONIAN 


W poprzednich częściach niniejszego opracowania przedstawione zostały ogólne 
teoretyczne podstawy opisu i analizy elementów układów optycznych w formali- 
zmie teorii grup Liego. Przedstawiono rozwiązanie zagadnienia propagacji przy- 
osiowego promienia świetlnego w światłowodzie gradientowym z dokładnością do 
wyrazów rzędu drugiego. Przedyskutowano także problem aberracji wywołanych 
uwzględnieniem w rozwinięciu hamiltonianu wyrazów czwartego rzędu. W tej 
części pracy rozważać będziemy przebieg promieni w światłowodzie gradientowym 
z dokładnością do wyrazów rzędu czwartego. Uwzględnienie tych wyrazów pro- 
wadzi do pojawienia się nieliniowości rzędu trzeciego w funkcjach opisujących bieg 
promieni. 

Zgodnie ze wzorem (I. 6) hamiltonian układu optycznego o współczynniku 
załamania n( p, q ,z) ma następującą postać: 


H(p,a,z)=-vn(p, q,2)?-p?. (1) 


*Praca finansowana z funduszy Problemu RR - I - 02 
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Załóżmy, że współczynnik załamania n zależy tylko od odległości od osi światłowodu: 
n = n(q?). Rozwinięcie wzoru (1) z dokładnością do wyrazów czwartego rzędu ma 
wtedy postać: 


1 p 2 1 y 
=> — pt +p + — 2 
ap? tug? + pt + 0g „żP 4 (2) 


Jak widać, hamiltonian składa się z dwóch części. Pierwsze dwa wyrazy opisują 
układ bezaberracyjny omówiony w części I. Pozostałe wyrazy — rzędu czwartego — 
są odpowiedzialne za występowanie aberracji. Ewolucję dowolnej funkcji zależnej 
od punktu w przestrzeni fazowej f( p, q ) opisuje zgodnie z (I. 20) następujący 
operator: 

exp(-zH). (3) 
Przebieg promienia świetlnego jest wyznaczony przez dwie funkcje: p (z)i q (2). 
Działanie operatora (3) na współrzędne p i q powoduje domieszanie” wyrażeń 
rzędu trzeciego, mamy więc do czynienia z problemem nieliniowym. Aby uwolnić się 
od nieliniowości, należy wprowadzić odpowiednie rozszerzenie przestrzeni fazowej. 
Ta nowa przestrzeń Z jest rozpięta na funkcjach połówkowych 1y1/2, 3x3/2 j 3y1/2 
zdefiniowanych przez wzory (II. 15, 16, 17). W ośmiowymiarowej przestrzeni Z 
operator ewolucji układu (3) działa liniowo, dzięki czemu obliczenia sprowadzają 
się do zwykłego mnożenia macierzy. Operator (3) jest szczególnym przypadkiem 
ogólnego elementu grupy A: G(v, w,u), przy odpowiednim wyborze parametrów 
v(z), w(z) i u(z). Po znalezieniu tych funkcji działanie operatora (3) na wektor 
fazowy (p , q ) obliczyć można ze wzoru (II. 19). 


2. REPREZENTACJA MACIERZOWA 


Hamiltonian (1) będący funkcją zmiennych fazowych p i q należy zastąpić 
operatorem A zdefiniowanym w równaniu (I. 7). Po wprowadzeniu operatorów 
"x7, związanych z funkcjami ” x?,, określonymi w (II. 1, 2, 3) operator H przyjmie 
postać: 


A 1 1 v 2v 
Ha 51 4 „gl RE PAZ A a UA 4 
Ino Xa ee gat ano? W ang Xo PX=2 + 3° Xo (4) 


Działanie tego operatora na punkt (*x1/2,3 x3/2,3 x1/2) rozszerzonej przestrzeni Z 
można przedstawić w postaci: 


1y1/2 ` A B B ]f}x!/2 
A 33/2 = 0 D 0 33/2 > (5) 
| Syt/2 | 0 0 8 A 
gdzie macierze A, B, C i D są następujące: 
4v ] 
0 2 3 0 
Be E O 735 p. Z 
no 3n2 
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> . 
0 0 4p 
nó 
0 


B = 1 u 0 1 (6) 
In} nd 
0 6v 0 
1° 
.— 0 4v 
no 2 
P> 0 -— 0 2v 
no 
0 0 TA 0 
no 


Macierzowym przedstawieniem operatora ewolucji (3) bedacego elementem grupy 


AS, która działa w przestrzeni funkcji określonych na rozszerzonej przestrzeni fazo- 
wej Z, jest macierz 


A B C 
exp(z| 0 D 6 la (7) 
0 0 A 


Jawna postać tej macierzy jako funkcji zmiennej z można znaleźć rozwiązując układ 
równań różniczkowych (II. 19), gdzie a(z) = (v(z), w(z), u(z)) jest zespołem para- 
metrów opisujących poszczególne rodzaje aberracji [1]. 

W niniejszej pracy macierz (7) obliczono bezpośrednio przez rozwinięcie funkcji 
wykładniczej i zsumowanie powstałego w ten sposób szeregu. W rachunkach wyko- 
rzystano wzór Perrona na n-tą potęgę macierzy. Dość żmudne obliczenia prowadzą 
w końcu do następującego wzoru: 


M~! M= y 2w M "' 


exp(—2H) = 0 D °/?(M `!) 0 ; (8) 
0 0 Mere? A 
gdzie: .5 s 
- K? 
meje AE 1-— KC, —S, 
7 ô —K? ng Sy i= KC; | 


1 3 1 
(0 +26 —=(nS4 + ez 
z ok 4 2) z(7 4 352) 


1 3 
-z ml — 452) — çz) el + 2C2) 


ab 
ER + 2C2) A (nfs + 482) — cz) | 
2yb 2b. 


1 » Y 
¿(MS + 352) zyj a ae 2C2) 


2v 
= —(——)z 
A (33) ; 
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a? 3a28 3a8? 83 
D m = a?y a(28y +a6) Plfy+2as) 836 
~ | ay? 6(206+ 48) 6(a6 + 218) Pe? |” 


7° 3976 378? 63 
gdzie: 
f2 
K= y = S = y2vno , 
no = 
_ v 1—2nop AV 6nop 
"= anż ya + rt =) 


A cola. 
kk kk 


Propagację promienia przedstawia trajektoria w zwykłej przestrzeni fazowej( p , q ). 
Ze skomplikowanego wzoru (8) wystarczy zachować tylko pierwszy wiersz, pozostałe 
bowiem nie niosą żadnej nowej informacji o przebiegu promienia: . 


PR 


= -1 (11/2 9,9/2 9,1/2 
ale). M (x Vy + 2w?y1/2) . 


3. ABERRACJE 


W części II pokazano, że każdy z sześciu parametrów (v, w) opisuje jakąś 
aberrację rozpatrywanego układu. Aberracje te są różne od aberracji znanych z 
analizy tradycyjnej, takich jak astygmatyzm, dystorsja, krzywizna pola itd., jed- 
nak wiążą się z nimi przez pewną transformację liniową. Ze wzoru (10) łącznie z 
wyrażeniami na macierze M `! i V , można odczytać następujące informacje: 


(a) w przybliżeniu rzędu pierwszego, tzn. bez uwzględnienia aberracji, ruch pro- 


een = : A 27 / 2no > 
mienla jest periodyczny z okresem rownym Ap = a TN F co wynika z 


postaci macierzy M 7?, w której występują funkcje sin(xz) i cos(xz); 
(b) występujące w układzie aberracje związane są z sześcioma parametrami vz, 
Vi, VO, U-1, V_2 Oraz w; ponieważ parametry te są wzajemnie niezależne, w 


układzie występować będzie w ogólności sześć różnych i niezależnych aberra- 
cji; 
(c) aberracje zależą od długości światłowodu, co wynika z postaci macierzy V , 


w której występują funkcje zależne od z; parametr w również zależy od z; 


*(d) niektóre aberracje mają część oscylacyjną z okresem Aa» = zAc bowiem 


wyrażają się przez funkcje sin(2kz), cos(2«z), sin(4«z) i cos(4xz); 


(1) 
[2] 


Bone € TOUHOCTLIO NO “YACHOB 4eTBEpToro nopzaxa. 
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(e) niektóre aberracje mają część zależną liniowo od z, współczynnik proporcjo- 
nalności jest równy ¢; 


(f) oscylacyjne części aberracji mogą zostać usunięte przez odpowiedni wybór 
długości światłowodu, która musi być wielokrotnością Agb; 


(g) przez odpowiedni dobór kształtu współczynnika załamania, to jest parametrów 
v ip, tak aby znikał współczynnik przy wyrazie liniowym (¢ = 0) można 
spowodować usunięcie aberracji, w których występuje część zależna liniowo 
od z; 


(h) jedyną aberracją, której nie można usunąć jest aberracja związana ze współ- 
czynnikiem w; w tradycyjnej nomenklaturze jest to pewna kombinacja krzy- 
wizny pola i astygmatyzmu. 


W przedstawionym wyżej formaliźmie teorii grup Liego możliwy jest jednolity opis 
powierzchni refrakcyjnych, a więc także powierzchni styku dwóch włókien o różnych 
parametrach. Możliwa jest także analiza światłowodów pozbawionych symetrii 
osiowej lub mających zmienną charakterystykę wzdłuż osi optycznej z. 
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SUMMARY 


Light ray propagation in a waveguide is analysed up to the fourth order. Six types of aberra- 
tion are discussed. 
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